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ANALYSIS OF THE ERRORS ASSOCIATED WITH THE DETERMINATION 


OF PLANETARY ATMOSPHERE STRUCTURE FROM MEASURED 

ACCELERATIONS OF AN ENTRY VEHICLE 

By Victor L. Peterson 
Ames Research Center 


SUMMARY 


An analytic method is developed for estimating the errors in the density 
and pressure structure of a planetary atmosphere constructed from measurements 
of accelerations experienced by an entry body* Among the sources of error 
considered both singly and collectively are the measuring accuracy of the 
accelerometers, the uncertainties in entry speedy entry angle, and aerodynamic 
coefficients, the frequency of data measurements, and the uncertainty of the 
attitude of entry bodies having nonzero lift-to-drag ratios* 

Results of sample calculations are presented to show the accuracy with 
which the extremes of a range of postulated model atmospheres for Mars can be 
defined from measurements of accelerations experienced by a spherically shaped 
entry body. It is found that if the accelerometer method were used only over 
that portion of the atmosphere traversed at speeds greater than sonic speed, 
then any of the atmospheres postulated for Mars can be defined reasonably well 
even when errors from all sources are combined in the most unfavorable manner. 
Some of the results obtained by the analytic method are compared with those 
from more precise numerical calculations and the agreement is found to be 
excellent . 


INTRODUCTION 


Properties of planetary atmospheres can be determined from observations 
on board an entry vehicle passing through the atmosphere* Experiments involv- 
ing a number of forms of observations, such as acceleration, spectral distri- 
bution and intensity of shock- layer radiation, oscillation frequency, etc., 
were proposed in reference 1 and the proposals were developed in more detail 
in reference 2. Considerable attention has been focused on one of these 
experiments in which knowledge of the acceleration experienced by a body, hav- 
ing known aerodynamic characteristics and entering an atmosphere on a ballis- 
tic trajectory, is used to construct the density and pressure structure of the 
atmosphere. A detailed procedure for obtaining atmosphere structure from 
measurements made with on-board accelerometers was given in reference 3j it 
was also shown that from the standpoint of flight mechanics, the ideal entry 
body for the atmosphere -structure experiment is an aerodynamically stabilized 
sphere. 



A number of sources of error in atmosphere definition are associated with 
the proposed method. Among these are the measuring accuracy of the acceler- 
ometers , the uncertainties in entry speed, entry angle, and aerodynamic coef- 
ficients ^ the frequency of data measurements, and the uncertainty in knowledge 
of the attitude of lifting entry bodies. Reference 2 gives an analytic method 
for predicting the possible inaccuracies in the definition of postulated Mars 1 
atmospheres due to just one of these sources, errors in measured accelerations. 
The method of reference 2 is limited, however, since it only allows for errors 
originating with the accelerometers and since assumptions of small errors and 
constant drag coefficients are inherent in the method. The effects of errors 
from this same source on the definition of the same postulated Mars atmos- 
pheres were estimated by a more precise numerical method in reference 3« 
Subsequent to these studies, general methods have been devised for assessing 
the effects of errors from the probable sources regardless of their magnitudes. 
These are presented herein and their use is demonstrated by some examples 
pertinent to the problem of determining the atmosphere of Mars. 


NOTATION 


A reference area for aerodynamic coefficients 

a^ resultant acceleration due to aerodynamic loads 

a s component of resultant acceleration due to aerodynamic loads directed 

along the flight path 

Cp coefficient of aerodynamic drag, drag/( l/2) pV 2 A 

Cp Q coefficient of aerodynamic drag evaluated at body attitude correspond- 

ing to zero lift 

Cp coefficient of aerodynamic lift, lift/( 1/2) pV 2 A 

Cp coefficient of resultant aerodynamic force, n/Cp 2 + c L 2 

g local acceleration due to gravity 

g Q acceleration due to gravity evaluated at surface of planet 

H n local atmosphere density scale height, defined by equation (4) 

h altitude above planet surface 

m mass of entry body 

p ambient pressure in atmosphere 

Re Reynolds number 

r distance from center of planet to mass center of entry body 
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t time 

V flight speed 

A( ) difference between quantity evaluated from measured data at a given 

instant of time and exact value of quantity evaluated at the same 
instant, (~)-( ) 

Q flight-path angle measured from local horizontal, positive down 

p ambient density in atmosphere 

a resultant angle of attack of entry body 

t dummy variable of integration 

9 tan -1 Cl/Cd 

(~) quantity evaluated from measured data 

Subscripts 

E quantity evaluated at entry 

F quantity evaluated at time corresponding to the final usable measure- 

ment of acceleration 

max maocimum value of function occurring in a specified entry 
min minimum value of function 


ANALYSIS 


Sources of error in constructing a planetary atmosphere from a time 
history of acceleration experienced by an entry body are identified and 
methods for assessing their influence on the desired results are developed 
in this section. Specifically, equations are developed for calculating the 
density, pressure, and density scale height variations with altitude. The 
equations are cast in forms that clearly show the roles played by the various 
factors in degrading precision of the results. 

The problem can be approached in a number of ways. For example, an 
error in the calculated density variation with altitude can be interpreted 
as being due entirely to errors in density at given altitudes, or entirely to 
errors in altitude at given densities. On the other hand, it might be inter- 
preted as the result of errors in density at given times in the trajectory 
plus errors in altitude at those times. The latter interpretation will be 
adopted here because it enables the prominent factors governing the problem to 
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be clearly identified. It also allows for the possibility of small errors in 
density and pressure at a given time, but yet for large errors in these quan- 
tities at a given altitude caused by a poor determination of altitude. It 
will be found that the equations for density, pressure, and altitude all have 
time as a common variable. The assumption is made that time is known exactly 
at the instants when acceleration is measured. The approach will be first to 
show how errors in density and pressure at a given instant of time are related 
to errors in knowledge of the trajectory at that same time. Then equations 
will be derived for the errors in trajectory variables in terms of time. 


Equations for Errors in Atmosphere Properties in Terms of Errors in 

Trajectory Variables 


An equation for the density p in terms of the drag coefficient Cp and 
the component along the trajectory of resultant acceleration due to aerody- 
namic loads can be written as' follows: 


P = 2 



a s 


( 1 ) 


Accelerometers located at the mass center of an entry body can be alined to 
measure components of acceleration along three mutually perpendicular axes. 

The vector sum of these components is the resultant acceleration due to aero- 
dynamic loads. Except when the lift of the entry body is zero, the magnitude 
and direction of this resultant vector are different from the component of 
acceleration along the flight path. This is illustrated in the sketch in 
which the longitudinal axis of the body is identified by the symbol x. The 

magnitude of the acceleration component 
along the flight path is given in terms 
of the resultant acceleration by the 
following equation 

a s = a^ cos cp 

where 

cos cp = ^ - — 

"Jl + (Cl/Cd) 2 

This result and equation (l) are combined to give the following expression for 
the density in terms of the resultant acceleration ap 

/ ^ \ a p cos cp , v 

9 = 2 (cy) w ^ 

To obtain the ambient pressure, the weight of a vertical column of atmosphere 
is integrated over the trajectory: 



f-t 

p = / gpV sin 6 dT 

Jo 


(3) 
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where it has been assumed that the altitude at time zero is high enough that 
the pressure at this instant is zero. The local density scale height is 
defined as 


H 



( 4 ) 


Quantities obtained from measured data, and therefore subject to error, 
will be denoted by a tilda over the symbol while those without a tilda will 
denote exact values. Equations for the ratios of measured values to exact 
values are formed by means of equations (2), (3), and ( 4 ) and are given below. 

£ 

P 


£d fR cosj)AA 2 
C D a R cos * W 


|j = - J gpv sin 0 dT / (5) 

Sfi = £ £ £ 

H p P P g 


A better form for these equations is obtained if a symbol is defined for the 
difference between values deduced from measurements and exact values. Let 


a( ) = n-( ) 


With this definition 



, equations 

1 + A cos cp 
cos cp 


1 

a v 


( 5 ) 


take the following form 




P 

P 


1 

P, 


gpV sin 0 



A cos cp^ ( - A sin e\ 

cos cp J sin 0 J 



dT 




^P = P P| 

H p p P s . 

J ( 6 ) 


Equations (6) relate the values of density, pressure, and scale height deduced 
from measurements to errors in acceleration, speed, path angle, gravity, lift- 
to-drag ratio, and coefficient of resultant force. The problem of determining 
the errors in trajectory variables, speed, path angle, and altitude from the 
measured resultant accelerations will be considered next. 
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Equations for Errors in Trajectory Variables 
The trajectory is governed by the following set of equations (ref. 3). 


dV 

dt 


+ 


v a e 

g dt 


dh 

dt 


^ (i pV 2 

m \2 


+ (Ii- 1 


-V sin 0 


- g sin 0=0 


cos 0 




J 


( 7 ) 


These equations written in terms of the resultant acceleration, a E , with the 
aid of equation (2) are 


a R cos cp - g sin 0 = 
+ f— - l) cos 0 + i 

w y g 

-V sin 0 


equations can be simplified if the flight-path angle is 
assumed to remain constant at the entry angle , that is, 0 = 0 E . The error 
resulting from this assumption diminishes with increasing entry angle and 
becomes zero for a $ 0 ° entry. It has been shown in references 1, 2, and 3 
that steep entries (50° to 90°) are preferred for this experiment so that 
assuming a constant path angle is not seriously restrictive, as verified by 
results of numerical calculations presented later. It is further assumed that 
the acceleration due to gravity is constant at the planet surface value. With 
these approximations, the equations for the trajectory deduced from measure- 
ments can be written as follows: 


Ct 

a R cos 9 


= 0 


( 8 ) 


<av + 

dt 

v de 

K dt 


dh 

dt 


The set of trajectory 


CLV Ai ^ 

+ a R cos cp - g Q sin 0 E = 0 


dh 

dt 


-V sin 0 E 




(9a) 


A similar pair of equations can be written for the trajectory which would be 
obtained under the same approximations (0 = 0 E , g = g Q ) if there were no 
errors in either a E or the entry conditions. The latter trajectory, which 
will be called the reference trajectory, is written as 
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I 


dV 

at 


ah 

at 


+ a R cos cp - g Q sin Qg = 0 

> 


= -V sin 0 e 


(9b) 


An equation for the difference AV between the computed speed and the 
reference speed is obtained when the first of equations (9b) is subtracted 
from the first of equations (9a)* This gives 


d(AV) 

dt 


+ % cos 9 - a-R - g Q A sin = 0 


Further algebraic manipulation puts the above equation into the following form 


+ (cos cp + A cos cp)AaR + A cos cp - g Q A sin Q-& = 0 (10) 

dt 

Equation (10) can be integrated analytically if cos cp, A cos cp, Aa^, and 
A sin are assumed constant for all time. The implications of making such 
assumptions will now be examined in detail. The assumption that cos cp and 
A cos cp are constant is exact for a spherical body for which Cp is always 
zero, for with Cp = 0, both cp and cp are identically zero. For nonspherical 
bodies undergoing planar oscillatory motion, the quantity cos cp will also be 
oscillatory with time, and for spinning nonspherical bodies undergoing coning 
motion, cos cp will not be oscillatory but rather will vary with the envelope 
of the total angle of attack. The present theory can be used to estimate the 
errors associated with the most extreme, assumptions for the cases involving 
nonspherical bodies. For example, if the rate at which the components of 
resultant acceleration were measured was not sufficient to allow the variation 
of angle of attack through the trajectory to be reconstructed, then the worst 
case, from the standpoint of errors, would be to assume that the angle of 
attack was identically zero for all time. This is equivalent to assuming that 
the lift-to-drag ratio associated with the measured data is zero so that 
A cos 9=1- cos cp. Now, for A cos cp to be constant, the quantity cos cp 
must be invariant. Again, the worst case can be assumed. For all bodies of 
revolution, except a sphere, the absolute value of the lift-to-drag ratio 
| Cl/Cd I increases from zero as the resultant angle of attack is either 
increased or decreased from zero. This resultant angle can never be larger at 
any point in the trajectory than the value at entry if the body is dynamically 
stable. The maximum possible value at entry would be known with reasonable 
precision so that, for any chosen probe configuration, an estimate could be 
made of the maximum value of |Cp/Ci)| (hence, 9) obtainable in the range of 
possible attitudes. A value so chosen gives a value of 9 denoted as 9max« 
Now consider the implication of assuming a constant error in measured acceler- 
ation Aap. Measurement accuracy of accelerometer systems is often given as 
a fraction of the maximum value measurable by the instrument, but it is 
likely that such errors would be random with time. For this analysis, however, 
they will be assumed constant at the maximum value. This means that the meas- 
ured accelerations would be biased from the true accelerations by a constant 
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maximum amount. This bias produces the largest errors in calculated atmos- 
pheres. Finally, assuming a constant error in path angle is consistent with 
assuming a constant path angle throughout the trajectory. The assumptions 
discussed above allow equation (10) to be rewritten in the following form: 

+ Aap> + (sec qJmax - l)a s - g Q A sin 0 E = 0 (ll) 

CLU 

Equation (ll) can be integrated analytically if the following substitution is 
made: 


a s = g 0 sin 0 E 


dV 

dt 


The result is 


AV = AV e + 


A sin 0 E -(sec c pm ax - l)sin 0 E 


Aa R 

So 


g 0 t - (sec cpmax - l)(V E 


- V) 

( 12 ) 


where 


sec cpmax 


" n/i + ( ^l/QD ) ]mx 


An equation for the difference Ah between the calculated altitude and 
the reference altitude at any instant of time is obtained when the second of 
equations (9b) is subtracted from the second of equations (9a). The result is 


d(Ah) 

dt 


-VA sin - AV(sin 0 e + A sin 0jj) 


(13) 


This equation can be integrated directly, with the aid of equation (12), to 
give 

A sin Oy . 

Ah = Ah F + ± (h - h F ) 

s in 


So 

2 


+ (sin 0 E + A sin 0 E ) 

^E “ ( sec 'ftnax ~ l)V E 


A sin 0 E - (sec c p^y - l)sin 0 E 


Aa R 
So . 


(t F - t) + 


( sec 'ftnax " l) 
sin 0 E 


(h - h F )j- 


(t F 2 - t 2 ) 
(14) 


Finally, equations (6) must be rewritten to be consistent with the 
assumptions just introduced* Two factors are affected* One is the ratio 
A cos cp/cos cp. This can be written 


A cos cp 
cos cp 


sec Tmax “ 


1 
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The other factor that requires explanation is [1 + (ACp/Cp)]. Since ffp must 
be taken equal to Cp 0 t o be consistent with the assumption that the body is 
in a nonlifting attitude. 


A . acd'N = c d 0 +ac d 0 

V C D J " C D 


(15) 


Since the aerodynamic characteristics of most simple bodies with thin shock 
layers are not strongly dependent on gas composition, it can be assumed that 
Cd q is a function only of Reynolds number and speed so that 

= CD Q (Re, V) - CDo(Re, V) 


Note that Cd in equation ( 15 ) should be evaluated at the body attitude 
selected for determining 9max» The equations for the density, pressure, and 
scale height become 


£ 

P 


£ 

P 


1 + 


r— j sec cpmax 
a R ' 


Cp 0 + ac Do ^ + avY 


Cd 

/ Aa R \ / A sin 0 R \ 

£° r%v Bln e E V - Y 3ec ^ 

P J o C Do + AC Do ^ ^ 

'or V v J 


V 


( 16 ) 


H 

_£ = 2 £ 

Hp p p 

J 

Equations (12), (l4), and (l6) constitute the central results of this 
analysis. They can be used to estimate the accuracy to which an atmosphere 
can be defined from measurements of resultant' acceleration. This is done as 
follows: First a model atmosphere is selected. Then, an accurate point-mass 

trajectory is computed with the selected model atmosphere and a given set of 
entry conditions. Finally, desired values for Aap, AVp, A sin Ahp, Tmaxj 
Cp, and Cp Q are selected and equations (l2), (l4), and (l 6) are used to cal- 
culate p, p, flp, and h. 


Note that total trajectory time tp plays an important role in this 
analysis so that the method of calculating the reference trajectory must be 
sufficiently precise to give reasonably accurate time histories of speed and 
altitude. Generally, the familiar analytic solution derived by Allen and 
Eggers (ref. 4) for steep entries with the assumption of zero gravity cannot be 
used. The reason is that in many situations the entry body will slow to ter- 
minal speed before impact with the planet surface and a considerable portion 
of the total trajectory time can be accumulated during this terminal descent. 
If gravity is neglected, the computed terminal speed can easily be in error by 
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more than a factor of 2, Under these circumstances, values of time to traverse 
given altitude changes would he in error. An accurate trajectory can be cal- 
culated either by solving equations (9b) for the reference trajectory (e.g., 
ref. 5) or by integrating the exact equations numerically. 


DISCUSSION 


Results given by the present technique will be illustrated by showing the 
accuracy with which extreme model atmospheres of the planet Mars can be 
deduced from acceleration measurements. The models of the Mars atmosphere 
given in appendix C of reference 6 will be used. To check their accuracy the 
approximate analytic results will be compared with results derived by the 
numerical method described in reference 3. The latter results are essentially 
exact, since no approximations are made either to the governing equations or 
their solutions. Only the density variation with altitude will be considered 
since it is the quantity least accurately determined from acceleration meas- 
urements. Since, from aerodynamic considerations, the optimum entry body 
shape for this experiment is an aerodynamically stabilized sphere (ref. 3), 
most of the discussion will pertain to this shape. In particular, the accu- 
racy with which atmosphere density as a function of altitude can be measured 
will be shown when errors in measured acceleration, entry speed, entry angle, 
and final altitude are assumed to occur both singly and combined in the most 
unfavorable fashion. Results will also be presented to indicate the effect of 

data sampling frequency on the descrip- 
tion of density variation with altitude. 
Finally, some estimates will be made of 
the penalties incurred when entry bod- 
ies other than spheres are used. 

Except where otherwise noted, a value 
of the ballistic parameter at entry 
(m/CpA)g of 0.25 slug/ft 2 will be 
assumed. Likewise, a typical entry 
speed of 26,000 ft/sec will be used. 

The coefficient of drag of the spheri- 
cal shape will be assumed to be inde- 
pendent of Reynolds number but will be 
allowed to vary with flight speed in 
the manner shown in figure 1 . The data 
given in figure 1 represent a mean 
fairing of sphere drag coefficient data 
for the range of Reynolds numbers 
appropriate to the entry conditions 
considered. 



Figure 1.- Assumed variation of the drag 
coefficient with flight speed for a 
spherical body. 


Results Obtained From Accelerations Measured Throughout the 

Entire Entry 

Consider first the effects of entry angle on the accuracy with which the 
extreme model atmospheres for Mars can be determined by means of accelerometers 
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with full-scale ranges chosen to match the specific entries. These are shown 
in figure 2 wherein biased errors in resultant acceleration equal to 0.1 per- 
cent of the maximum values experienced have been assumed. It is shown for 



Figure 2.- Effect of entry path angle on the definition of extreme model atmosphere for Mars: 
spherical body; (m/Cj>A) E =0.25 slug/ft 2 ; V R = 26,000 ft /sec; biased acceleration error 
Aa R = -0.001 (aR^max* 

these circumstances that entry angle affects the accuracy of atmosphere 
determination only slightly but that results are quite different, depending 
upon the model atmosphere encountered. The accuracy with which the maximum 
scale height atmosphere can be determined is poorer than for the minimum 
atmosphere, primarily because the flight time through the former is so much 
longer than through the latter. Errors are accumulated over a time period 
nearly six times longer in the maximum atmosphere than in the minimum atmos- 
phere. Note that the method gives results in excellent agreement with those 
obtained numerically from the exact equations. 

The results shown in figure 2 are idealized in the sense that the accel- 
erometer ranges are matched to the chosen conditions while, in practice, it 
may not be known which atmosphere within the possible extremes will be 
encountered on a first mission into the Martian atmosphere. This means that 
the accelerometers will have to be designed to accommodate any magnitude of 
acceleration that can possibly be experienced. The maximum possible accelera- 
tion will occur during entry into the minimum scale height atmosphere with the 
steepest flight-path angle. It has a value of about 438 times the gravita- 
tional acceleration at the surface of Mars. For illustrative purposes, it 
will be assumed that the accelerometers must be designed to measure 
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accelerations as large as 450 g Q with an accuracy of 0.1 percent of this 
value (0.45 go)* The most unfavorable situation from the standpoint of atmos- 
phere definition clearly would be that of designing for the maximum possible 
accelerations and actually experiencing the minimum values . The minimum 
values correspond to entering the maximum scale height atmosphere with the 
shallowest path angle. Therefore, emphasis will be placed on these conditions. 

Figure 3 shows the additional penalty of designing for the largest accel- 
erations in the minimum atmosphere and actually encountering the maximum 
atmosphere. The additional penalty in atmosphere definition is unacceptable. 

In order to gain insight into what steps might be taken to improve matters, 
this extreme case will be examined in more detail. This can be done with the 
aid of the results of figure 4. These results show how the error in atmos- 
phere definition is divided among errors in measured accelerations, errors in 
speed due to integrating these acceleration errors, and errors in altitude due 
to integrating the errors in speed. It is shown that the atmosphere can be 
described quite well if speed and altitude are measured accurately by an inde- 
pendent method. The definition is equally good over a less extensive range of 
altitude if only the altitude is determined by an independent method. The 
altitude at which speed errors begin to increase rapidly corresponds to a 
speed of about 2000 ft/sec. The largest errors are shown to result from 
errors in altitude. Although the altitude errors are large over the entire 
trajectory, they are, in fact, accumulated largely over the low-speed part of 
the trajectory where speed errors are large. The high-altitude results are 
degraded by these errors only because the altitude at entry is not known a 
priori so that changes in altitude are, in effect, determined relative to the 
final altitude which is assumed to be zero at probe impact. These considera- 
tions suggest that either the accelerometer accuracy must be improved consid- 
erably or else different methods for measuring atmosphere properties must be 
used once the probe speed is reduced to low values. The improvements in 
atmosphere definition that can be expected to result from both of these tech- 
niques will be discussed next. 

The effects of increasing accelerometer accuracy by factors of 2 and 10 
are shown in figure 5* and the effects of terminating the accelerometer exper- 
iment at speeds ranging from that for impact at 378 ft/sec up to 1000 ft/sec 
are shown in figure 6. The results in figure 5 show that about an order of 
magnitude increase in accelerometer accuracy (to 0.01 percent of the maximum 
measurable) would be needed to provide good definition of the atmosphere at 
all altitudes. This accuracy is believed to be beyond that obtainable with 
either single or dual range accelerometer systems. The most attractive prac- 
tical means for improving the definition of the atmosphere at altitudes where 
accelerations are large appears to be to accept data only over that range of 
altitudes for which the speed is greater than sonic speed. The results of 
reference 2 indicate that speed, altitude, and atmosphere properties can be 
determined well by direct measurement at subsonic speeds. For the remainder 
of this discussion, it will be assumed that the accelerometer experiment would 
be terminated at a speed of 750 ft/sec. Of course, errors in the quantities 
measured at low speed by the methods described in reference 2 will affect the 
accelerometer experiment because of its dependence on altitude at the final 
time - in this case, that corresponding to a flight speed of 750 ft/sec. 
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igure 6.- Effect of speed at which accelerometer experiment is terminated or 
a model Mars atmosphere; spherical body; (m/Ci)A) E = 0.25 slug/ft 2 ; V E = 2 
« 50°: biased acceleration error Aap = -0.45 g 0 . 



Effects of errors in the low speed measurements can be included in this study 
by the assumption of various errors in the final altitude. 


Results Obtained From Accelerations Measured Only at Speeds 
Greater Than 750 Ft/Sec 

First, attention will be directed at an examination of the individual and 
combined effects on atmosphere definition caused by less stringent accelerom- 
eter accuracies and errors in entry speed, entry angle, and altitude at the 
final speed. Then, consideration will be given to the effects of data sam- 
pling frequency. Finally, estimates will be made of the errors caused by the 
use of nonspherical entry bodies. 

Effects of errors in measured acceleration, entry speed, entry angle, and 
final altitude .- The effects of reducing accelerometer accuracy by factors of 
2 and 3 are illustrated in figure 7 . Over most of the altitude range the 
effects vary linearly with the change in accuracy and are about the same mag- 
nitude for both positively and negatively biased errors. The comparison with 
results from numerical calculations is excellent throughout. The effects of 
errors in entry speed, entry angle, and final altitude are shown in figures 8 , 
9, and 10, respectively. In each case it was assumed that accelerations were 
known exactly, that is, AaR =0. It is seen that degradation of atmosphere 
definition resulting from reasonable values for each of these sources is not 
overwhelming, and in each case, the analytic method provides excellent predic- 
tions. It is noted that results of figure 10 show errors in altitude through- 
out the trajectory to be equal to the errors assumed to exist at the final 
speed. This means that the inaccuracy of the accelerometer experiment due to 
errors in various methods for measuring the atmosphere directly at low speed 
can be readily determined since altitude is the common link between the high 
and low speed methods . 

The quality of atmosphere definition obtained when positive and negative 
errors in measured accelerations, entry speed, entry angle, and final altitude 
are combined in the most unfavorable manner is shown in figure 11. As 
expected, if the minimum scale height atmosphere were encountered, it would be 
defined more accurately than would the maximum scale height atmosphere should 
it be encountered instead. In either case, the present uncertainty of the 
Mars atmosphere would be reduced considerably by the measurements. Results of 
calculations not presented herein showed the data of figure 11, which are for 
values of speed and m/CpA at entry of 26,000 ft/sec and 0.25 slug/ft 2 , 
respectively, to be essentially the same for values of entry speed ranging 
from 22,000 to 28,000 ft/sec and values of m/CpA at entry ranging from 0.20 
to O .35 slug/ft 2 , provided the ranges for the accelerometers were adjusted 
for each combination of Ve and m/CpA. 

Effect of data sampling frequency .- The effects of data sampling 
frequency cannot be easily evaluated by analytic techniques. This possible 
source of degradation of the accelerometer measurements arises from the prac- 
tical necessity of minimizing the amount of data transmitted from a probe to 
Earth. All of the results presented so far were based on the assumption that 
the acceleration measurements were described by a sufficient number of points 
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Figure 9.- Effect of errors in entry angle on the d< 
spherical body; (m/C;QA) E = 0.25 slug/ft 2 ; V E = 2 ( 
V F = 750 ft/sec. 

700- 
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Figure 10.- Effect of errors in final altitude on the 
spherical body; (m/Cr>A) E = 0.25 slug/ft 2 ; V E = 26,1 
V F « 750 ft/sec. 





Figure 12.- Effect of time interval between measured accelerations on the definition of extreme 
model atmosphere for Mars: spherical body (m/QDA) E = 0.25 slug/ft 2 ; Vg = 26,000 ft/sec; 

(VjjO^ = 750 ft/sec; biased acceleration error Aa^ = -0.45 g 0 j 0 E = 50° and 90° for maximum 
and minimum atmospheres, respectively. 



to be effectively infinite. In other words, increasing the amount of data 

could not improve the accuracy. The effect of reducing the frequency of data 

sampling to allow for as much as 2 seconds of time between data points is 
shown in figure 12 . These results were obtained by means of the numerical 
method of reference 3. Second-order interpolation on the logarithm of meas- 
ured acceleration as a function of time was used. Note that these calcula- 
tions also include a representative combination of errors from other sources. 
The magnitude of the sampling intervals considered has very little effect on 
the quality of atmosphere definition except at the very low density levels 

where the atmosphere is first sensed by the instruments and at very low alti- 

tudes in the minimum atmosphere case. These results are considered remarkable 
in view of the fact that the acceleration time history is defined by only 12 
data points for the 90° entry into the minimum scale height atmosphere when 
the data sampling frequency is 2 seconds. 
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Figure 13.- Aerodynamic characteristics 

obtained from Newtonian impact theory for 
various entry body shapes. 


Effects of using entry bodies 
other than a sphere .- Estimates have 
been made of the errors in atmosphere 
definition caused by the use of entry 
bodies other than spheres. These esti- 
mates were made for the shapes illus- 
trated in figure 13. The aerodynamic 
characteristics presented in figure 13 
were obtained from Newtonian impact 
theory. It is obvious from the data in 
figure 13 that conclusions regarding 
the relative merits of the various 
entry body shapes will depend strongly 
on the angle of attack selected for the 
comparison. In this discussion, a com- 
mon angle of attack of 20° has been 
chosen; that is, the drag coefficient 
and the lift-to-drag ratio parameter 
cpmax were evaluated at this angle of 
attack for each configuration. These 
values, were used in equations (12), 
(l 4 ), and (l6) to obtain the results 
presented in figure l 4 . The aerody- 
namic characteristics were assumed 
invariant with Reynolds number and 
speed. No errors from any sources 
other than that related to unknown 
angle of attack were included in these 
calculations. The results of figure l 4 
show that the penalty in atmosphere 
definition associated with the angle of 
attack uncertainty makes two of the 
bodies definitely not suited for this 
application; however, the penalty for 
the body most closely resembling a 
sphere is not excessive. It will be 
recalled that all of the assumptions 
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made in the analysis, for deriving these results, were such as to make the 
effects of errors extreme. Actually, the errors due to using nonspherical 
configurations should not be as large as those shown in figure l4 for the same 
entry conditions. Also, the possibility exists of improving accuracy by cor- 
recting measurements for angle of attack. The body attitude can, in principle, 
be inferred from three -component accelerometer data and known aerodynamic 
characteristics. This possibility is discussed in more detail in references 2 
and 3 . 


CONCLUDING REMARKS 


An approximate method has been developed for estimating the errors in 
the density and pressure structure of a planetary atmosphere constructed from 
measurements of accelerations experienced by an entry vehicle. Among the 
sources of error considered both singly and in combination are accelerometer 
inaccuracies, uncertainties in entry speed, entry angle, final speed and 
altitude, aerodynamic characteristics, and frequency of data measurements. 

The effects of uncertainty of the attitude of entry bodies having nonzero 
lift -to- drag ratios were also considered. 
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The method was applied to the problem of estimating the accuracy with 
which the extremes of a range of postulated atmospheres for Mars could be 
defined by measurements of the accelerations experienced by a spherically 
shaped entry body. It was concluded that use of the accelerometer method for 
obtaining atmosphere structure should probably be confined to the portion of 
the atmosphere traversed at speeds greater than sonic speed. For this situa- 
tion, it was shown that any of the atmospheres within the extremes postulated 
for Mars can be defined reasonably well even when errors from all sources are 
assumed to be combined in the most unfavorable manner. Selected portions of 
the results were compared with results obtained from more precise numerical 
procedures and the agreement was found to be excellent# 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., April 13, 1965 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of ^phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate disseminatio?t 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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